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Field of the Invention 

The present invention concerns methods and apparatus for the identification 
and characterization of regions of different mechanical properties in a target media, 
25 such as a biological tissue. 



Background of the Invention 

The early detection of diseases, such as cancer and heart disease, may 
30 significantly improve patient survival. For example, present methods of breast cancer 
detection include screening mammography and palpation, either by patient self- 
examination or clinical breast exam. Palpation relies on the manual detection of 
differences in tissue stiffness between breast lesions and normal breast tissue. The 
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success of palpation is due to the fact that the elastic modulus (or Young's modulus) 
of malignant tumors is often an order of magnitude greater than that of normal breast 
tissue. That is, cancerous lesions feel "hard" or "stiff 5 as compared to normal breast 
tissue. See T. Krouskop et al., Ultrasonic Imaging 20, 260-274 (1998); A. 
5 Saravazyan et al., Acoustic Imaging 21, 223-240 (1995). 

As another example, atherosclerosis is a medical condition that involves the 
stiffening of the arteries, which can occur prior to and in addition to the formation of 
focal lesions. A human artery has three layers. The innermost layer is relatively thin 
compared to the other two layers and is commonly referred to as the endothelium 

10 layer or intimal layer. The middle layer is a smooth muscle layer (media layer), and 
the outer layer is a connective tissue layer (adventita). Atherosclerosis can be caused 
by plaque build up between layers in the arteries, which can lead to conditions that 
include blockages in blood flow, poor circulation, myocardial infarction, aneurisms, 
and stroke. The development of atherosclerosis can be influenced by factors such as 

15 diet, exercise, smoking habits, and other medical conditions such as diabetes. 
Treatments for atherosclerosis include drug therapy, changes in diet, increased 
exercise programs, and smoking cessation. 

Focal lesions, which protrude into the blood vessel, may be detected by a 
variety of methods, including cardiac catheterization, intravascular ultrasound, 

20 conventional B-mode and Color Doppler ultrasound, and electron beam computed 
tomography. However, atherosclerosis may occur prior to the formation of focal 
lesions and does not necessarily involve plaque protrusions that could be detectable 
using conventional testing methods. Atherosclerosis may also be characterized by 
increased blood pressure. However, blood pressure levels do not directly correlate to 

25 degrees of atherosclerosis, and elevated blood pressure levels may not exist in a 
subject without atherosclerosis. 

U.S. Patent No. 5,921,928 to Greenleaf uses sound waves to vibrate tissue and 
monitors the response of the vibrating tissue. A disadvantage of such an approach is 
that different transmit and detection systems may be required, and multiple pushing 

30 cycles at a single location may be required. Moreover, ultrasonic stimulated acoustic 
emissions may not result in a form feasible for actual clinical diagnosis. 

Summary of the Invention 
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According to some embodiments of the invention, methods, systems, and 
computer program products are provided. In some embodiments, a pushing pulse is 
delivered from an ultrasound transducer array having a plurality of elements to a 
target region within a medium to displace the target region to a displaced position. A 
5 tracking pulse is delivered from the ultrasound transducer array to the target region, 
and a plurality of tracking signals are received from locations in the target region, 
each tracking signal being responsive to the tracking pulse. The received tracking 
signals can be received in parallel receive mode, including parallel receive mode in 
the time domain and/or frequency domain, and can include parallel beam formed 
10 signals. 

In further embodiments, a first tracking pulse is delivered from an ultrasound 
transducer array to a target region within a medium. A first set of tracking signals 
from locations in the target region and responsive to the tracking pulse in the target 
region is received at the ultrasound transducer array to detect an initial position for the 

15 target region. A pushing pulse is delivered from the ultrasound transducer array to the 
target region to displace the target region to a displaced position. A second tracking 
pulse is delivered from the ultrasound transducer array to the target region. A second 
set of tracking signals from locations in the target region is received responsive to the 
second tracking pulse at the ultrasound transducer array to detect the displaced 

20 position of the target region. Delivering a first tracking pulse, receiving a first 

tracking signal, delivering a pushing pulse, delivering a second tracking pulse, and 
receiving a second tracking signal can be repeated sequentially to provide a series of 
cycles. The pushing pulses can be delivered to different target regions during the 
series of cycles to provide a plurality of displaced positions. 

25 In some embodiments a blood vessel and/or cardiac tissue in a subject can be 

evaluated by detecting a mechanical property of a vessel wall and/or cardiac tissue to 
provide at least first and second values associated with the mechanical property at at 
least two timing points during a cardiac cycle of the subject. The first and second 
values can be compared, for example, to evaluate the cardiac health of the subject. 

30 In further embodiments, a blood vessel and/or cardiac tissue in a subject can 

be evaluated by detecting a first value associated with a mechanical property of a 
vessel wall and/or cardiac tissue at a first spatial point within a first layer of the vessel 
wall and/or cardiac tissue. A second value associated with the mechanical property at 
a second spatial point of the vessel wall and/or cardiac tissue within a second layer of 
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the vessel wall and/or cardiac tissue can be detected. The first and second values can 
be compared. 

While the invention has been described above primarily with respect to the 
various method aspects of the invention, both systems and/or computer program 
5 products are also provided. 

Brief Description of the Drawings 
Figure 1 is a schematic diagram of systems for evaluating blood vessels using 
an ultrasound transducer array and parallel receive mode processing according to 
10 embodiments of the present invention. 

Figure 2 is a schematic diagram of systems for evaluating a target region 
using ultrasound interrogation techniques according to further embodiments of the 
present invention. 

Figure 3 is a block diagram illustrating further systems using ultrasound 
15 interrogation techniques according to further embodiments of the present invention. 

Figure 4 illustrates the sequential scanning of the two dimensional plane of 
target regions of Figure 3, with different forced regions in each cycle, to produce a 
two dimensional displacement map for each forced region. 

Figure 5 illustrates the signal processing operations implemented by the signal 
20 processing device (31) of Figure 3, in which the two dimensional displacement maps 
for each forced region, as generated in Figure 4, are combined into a single image. 

Figure 6 A, is an example remote palpation image of a breast lesion phantom. 
Figure 6B is the corresponding B-mode image of the breast lesion phantom 
shown in Figure 6A. 

25 Figure 7 illustrates embodiments of an orientation of a transducer and the 

associated different dimensions: axial, azimuthal, and elevation. 

Figure 8 illustrates embodiments of generating a three-dimensional volume 
using a translation stage connected to the transducer, which can allow the 
interrogation of multiple axial/azimuthal planes by translating the transducer in the 
30 elevation dimension. 

Figure 9 illustrates embodiments of a two-dimensional transducer (which has 
several rows of elements) to interrogate a three-dimensional volume. 
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Figure 10 is a schematic diagram of systems using parallel receive mode 
processing according to embodiments of the present invention. 

Figure 11 is a flow chart illustrating operations according to embodiments of 
the present invention. 

5 Figure 12 is a schematic diagram of systems for interrogating blood vessels 

and cardiac tissue according to embodiments of the present invention. 

Figures 13-16 are flow charts illustrating various operations according to 
embodiments of the present invention. 

Figures 17 and 18 are images obtained using systems according to 
10 embodiments of the invention; 

Figure 19 is a displacement value plot obtained from the ultrasound 
measurements of Figures 17 and 18 according to embodiments of the present 
invention. 

15 Detailed Description of the Preferred Embodiments 

The present invention now will be described more fully hereinafter with 
reference to the accompanying drawings, in which illustrative embodiments of the 
invention are shown. This invention, however, should not be construed as limited to 
the embodiments set forth herein; rather, these embodiments are provided so that this 

20 disclosure will be thorough and complete, and will fully convey the scope of the 
invention to those skilled in the art. 

Embodiments of the present invention can use interrogation techniques to 
characterize a target region, such as to identify regions of greater or lesser stiffness in 
biological tissue. Ultrasound interrogation techniques used to interrogate blood 

25 vessels using receive mode parallel signal processing are first described. Then, 
various interrogation techniques applied to a target region, such as biological tissue, 
are described. Then, receive mode parallel processing techniques aire described in 
greater detail. Various techniques used to study blood vessels are discussed. 

30 Receive mode parallel processing ultrasound techniques to interrogate blood vessels 

According to embodiments of the present invention described with reference 
to Figure 1, ultrasound techniques using receive mode parallel processing can be used 
to interrogate blood vessels. It will be understood that various elements and 
techniques can be used independently of other elements and techniques illustrated in 
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Figure 1. For example, ultrasound interrogation techniques can be used with or 
without receive mode parallel processing. Ultrasound interrogation techniques can 
also be used to interrogate other types of tissue, such as cardiac tissue, muscle tissue, 
kidneys, breast tissue, or other biological tissue. 
5 Referring to Figure 1, an ultrasound array 1010 for interrogating a target 

region 1016 of a tissue medium 1012 that includes a blood vessel 1004 is provided. 
The ultrasound array 1010 includes transducer elements 1010A-1010H connected to a 
controller 1000. The controller 1000 can control pulses produced by the transducer 
elements 1010A-1010H and can include a signal analyzer 1002 for analyzing 

10 ultrasound signals received by the transducer elements 1010A-1010H. The signal 
analyzer 1002 can use receive mode parallel processing techniques to analyze signals 
from array elements 1010A-1010H. Any number of array elements 1010A-1010H 
can be provided, including a single transducer and one or two dimensional transducer 
array configurations. It should be understood that embodiments of the present 

15 invention are not limited to the configuration shown in Figure 1. Certain elements 
shown in Figure 1 can be omitted or the functionality of certain elements can be 
combined with other elements. For example, the signal analyzer 1002 can be 
provided as part of the controller 1000 or the signal analyzer 1002 can be provided as 
a separate component from the controller 1000. 

20 As shown in Figure 1, each element 1010A-1010H in the array 1010 can 

produce an ultrasound pulse 1014 that propagates through the tissue medium 1012 to 
the blood vessel 1004 in the target region 1016. The ultrasound pulse 1014 creates a 
reflected ultrasound wave. In some instances, the ultrasound pulse 1014 can be 
delivered at an intensity sufficient to cause movement of the target region 1016. 

25 When the ultrasound pulse 1014 reaches a reflective boundary between tissues (e.g. 
between fluid and soft tissue, soft tissue and bone, etc.), some portion of the sound 
waves are reflected back to the ultrasound array 1010, whereas a remaining portion of 
the pulse propagates potentially reaching another reflective boundary. Ultrasound 
waves that may be reflected from within the tissue medium 1012 can be received as a 

30 signal by the ultrasound array 1010. As shown in Figure 1, the reflected waves can 
be received by one or more array elements 1010A-1010H along exemplary directions 
1018. Accordingly, the reflected signals can be received at the array elements 1010A- 
101 OH in a plurality of positions, such as shown, for example, with respect to array 
elements 1010D and 101 OH. The signals can be received from locations within the 
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target region 1016. The signal analyzer 1002 analyzes the signals received by the 
ultrasound array 1010, for example, to calculate the distances from the array 1010 to 
the point at which the ultrasound waves were reflected and the intensities of the 
various reflected waves. 
5 As illustrated for clarity and ease of representation, the pulse 1014 is produced 

by element 101 0B and the reflected ultrasound wave directions 1018 illustrate signals 
received by array elements 1010D and 101 OH. However, it should be understood that 
each element 1010A-1010H can produce a propagation pulse, such as pulse 1014, 
and/or receive reflected wave signals, such as waves propagating along directions 

10 1018. Accordingly, the pulse 1014 produced by element 1010B may result in 
reflected wave signals that can be received by any combination of element(s) 1010A- 
1010H. Moreover, in some embodiments, a single transducer element can produce a 
pulse and, in turn, receive the reflection wave signal. 

As shown in Figure 1, the pulse 1014 can be a de-focused beam with respect 

15 to the target region 1016. The reflected waves can be received from various locations 
in the target region 1016 and processed in parallel, for example, using a time-delayed 
signal to take into account varying distances between the origin of the reflected wave 
and each ultrasound transducer element 1010A-1010H. The time delay can be 
configured so that the reflected waves recorded by each transducer are in phase before 

20 they are summed and/or analyzed. Thus, the resulting signal can be a sum of focal 
delayed signals received at each of the transducer elements 1010A-1010H. Multiple 
sums of differently focused delayed signals can be used for a single pulse. 
Accordingly, increased information may be obtained for a single pulse, which can 
result in increased frame-rates and/or reduced tissue heating. 

25 In some embodiments, a pulse or a series of pulses can be produced by various 

elements 1010A-1010H in the ultrasound array 1010 to remotely characterize tissue 
stiffness, such as by using a pulse with sufficient ultrasonic radiation force so as to 
displace the target region 1016. One particular approach is referred to as remote 
palpation. In remote palpation, acoustic radiation force is used to apply localized 

30 forces within tissue, and the resulting tissue displacements are mapped using either 
ultrasonic correlation based methods, or other pattern matching methods. A volume 
of tissue that is stiffer than the surrounding medium (i.e., a lesion) may distribute the 
force throughout the tissue beneath it, resulting in larger regions of displacement, and 
smaller maximum displacements. Remote palpation is described in K. Nightingale, 
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Ultrasonic Generation and Detection of Acoustic Streaming to Differentiate Between 
Fluid-Filled and Solid Lesions in the Breast (Ph.D. Thesis, Duke University, 
September 1997), and in K Nightingale et al., Proceedings of the 1997 IEEE 
Ultrasonics Symposium, 1419-1422 (Toronto,. CA, Oct 1997). Another approach is 
5 referred to as acoustic radiation force impulse ("ARFI"), which uses ultrasound pulses 
to push on tissue, and subsequently measures the mechanical response of the tissue. 
Various ARFI techniques are disclosed in U.S. Patent No. 6,371,912 to Nightingale et 
al., the disclosure of which is hereby incorporated by reference in its entirety. 

In some embodiments of the present invention using ARFI techniques, a 

10 "pushing" pulse can be delivered from the ultrasound array 1010 to displace the target 
region 1016 such that the displacement of the target region can be detected. A 
"tracking" pulse can be delivered from the ultrasound transducer array 1010 to the 
target region 1016. The tracking pulse and/or the pushing pulse may be represented 
by pulse 1014 as shown in Figure 1. A tracking signal, such as the reflected 

15 ultrasound wave, can be received by the ultrasound array 1010 from a plurality of 
directions 1018. If a tracking pulse is broadly focused, a plurality of tracking signals 
can be received at the transducer array 1010. The tracking signal can be used to 
detect the displaced position of the target region 1016. Regions of altered stiffness 
can be detected, for example, based on characteristics of the displacement at various 

20 locations detected in the target region 1016. In some embodiments according to the 
invention, the tracking signals received from locations in the target region 1016 and 
responsive to the tracking pulse may be processed using receive mode parallel 
processing techniques. For example, the tracking signals can be parallel beam formed 
signals. Time domain and frequency domain receive mode parallel processing 

25 techniques can be used. For example, multiple sums of differently focused delayed 
signals received at each of the plurality of elements in the ultrasound transducer array 
can be used. Accordingly, receive mode parallel processing techniques can result in 
increased frame-rates and/or reduced tissue heating. 

In some embodiments, an initial tracking pulse can be delivered from the 

30 ultrasound transducer array 1010 to the target region 1016 prior to the pushing pulse. 
An initial tracking signal responsive to the initial tracking pulse can be received to 
detect an initial position for the target region. A difference between the initial 
position and a displaced position can be used to detect regions of varying degrees of 
stiffness. For example, a region of decreased displacement relative to other regions 
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can indicate that the region of decreased displacement has a greater stiffness relative 
to other regions. A region of lesser stiffness can be detected by an indication that the 
region has greater displacement relative to other regions. 

The tracking pulse(s) and the pushing pulse can be delivered at different 
5 intensity levels so that the intensity level of the pushing pulse is greater than the 
intensity level of the tracking pulse. For example, the intensity level of the pushing 
pulse can be between about 1.0 or 10.0 W/cm 2 and about 1000.0 or 10 5 000.0 W/cm 2 , 
and the intensity level of the tracking pulse can be less than about 1.0 W/cm 2 or about 
0.72 W/cm 2 . The pushing pulse can be delivered for between about 0.025 or lower to 

10 about 0.5 milliseconds or more (for example, up to about 10 milliseconds) at various 
intensity levels. The intensity level can be selected based on the length of the 
delivered pulse. In some embodiments, transmitting a tracking pulse and receiving a 
tracking signal can be repeated for a single pushing pulse. For example, an initial 
position can be detected by an initial tracking pulse prior to the pushing pulse, and a 

15 plurality of displaced positions can be detected by a plurality of subsequent tracking 
pulses after the delivery of the pushing pulse. 

As shown in Figure 1, the target region 1016 includes a blood vessel 1004. 
Other media can be used, including biological tissue such as breast tissue, muscle 
tissue or cardiac tissue. Embodiments according to the present invention can be used 

20 to detect a region of greater stiffness in biological tissue including tumors, hardened 
blood vessels, or muscle tissue with greater muscle tone. 

For example, blood vessels (such as blood vessel 1004) and/or cardiac tissue 
in a subject can be evaluated by detecting a value associated with a mechanical 
property of a vessel wall and/or cardiac tissue to provide at least two values associated 

25 with the mechanical property at two or more timing points during a cardiac cycle of 
the subject. The values associated with the mechanical property can be compared, for 
example, to evaluate cardiac and/or vascular health. As another example, values 
associated with a mechanical property of a vessel wall and/or cardiac tissue can be 
detected at spatial points in different layers in a vessel wall and/or cardiac tissue. The 

30 values of the mechanical property at the different spatial points can be compared to 
evaluate cardiac and/or vascular health. The values associated with the mechanical 
properties can be detected using ultrasound techniques as described herein. 
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Ultrasound interrogation techniques 

Further embodiments according to the present invention are shown in Figure 
2. A controller 1025 including a signal analyzer 1021 is connected to an ultrasound 
transducer 1020. The controller 1025 can control the transducer 1020 to deliver a 
5 pulse along a propagation path 1024 to a target region 1026 in a tissue medium 1022. 
A reflected wave can travel along the propagation path 1024 to be received by the 
transducer. The pulses can include "tracking" or "pushing" pulses, and the reflected 
wave can be a tracking signal. 

As shown in Figure 3, an exemplary ultrasound transducer array 20 is 

10 provided. The ultrasound transducer array 20 may be a one-dimensional array or a 
two-dimensional array. The ultrasound transducer array 20 may be contacted to a 
target medium, such as a tissue medium 10. The array can be directed to a two- 
dimensional plane comprising one or more, and in some embodiments, a plurality of, 
target regions 11 within the tissue medium. A transmit circuit 21 may be operatively 

15 associated with the transducer array to deliver high energy "pushing" pulses to a 
forcing region among the target regions (i.e., pulses that can induce a physical 
displacement of the tissue within the target regions), as well as for delivering 
relatively lower energy "tracking" pulses. A receive circuit 22 can be connected to 
the transducer array 20 to receive information from the target regions 11 for 

20 subsequent signal processing. The transmit circuit 21 and detector circuit 22 may 
both be operatively associated with an appropriate control circuit 23 that triggers the 
pushing pulses and tracking pulses, organizing information received from the target 
regions for subsequent signal processing, and which also cycles the pushing pulses 
and corresponding tracking pulses through different forcing regions. 

25 Information received by receive circuit 22 can be stored in a memory device 

30 such as a random access memory or other suitable memory device, which serves as 
a medium for storing both initial and displaced positions of target regions. A signal 
processing device or signal processor 31 is operatively associated with the memory 
device 30, and can generate initial images for particular forced regions and a single 

30 combined image for a plurality of forced regions. 

According to embodiments of the present invention, one or more of the 
following functionalities may be provided: 

(a) delivering a set of tracking pulses from a plurality of transducer elements 
in an ultrasound transducer array to one or a plurality of target regions in a two- 
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dimensional plane within said medium to detect an initial positions for said one or 
plurality of target regions; 

(b) storing said initial positions for said one or plurality of target regions; then 

(c) delivering a first set of pushing pulses from said plurality of transducer 
5 elements to a forcing region among said target regions to displace said target regions 

to subsequent (e.g., displaced) positions; 

(d) delivering a second set of tracking pulses from said plurality of transducer 
elements in said ultrasound transducer array to said one or plurality of target regions 
to detect subsequent positions for said one or plurality of target regions, 

10 (e) storing said displaced positions for said one or plurality of target regions; 

(f) repeating steps (a) through (e) in a series of cycles, with said pushing and 

tracking pulses being delivered from a different plurality of transducer elements or the 

same plurality of transducer elements in said array to a different forcing region, and 

preferably to a plurality of different target regions, during each of said cycles; 
15 (g) generating a two-dimensional displacement map from each of said initial 

positions and displaced positions for each of said forcing regions to produce a 

plurality of two-dimensional displacement maps; and then 

(h) combining said plurality of two-dimensional displacement maps into a 

single combined image, with a region of increased stiffness being indicated by a 
20 region of decreased displacement within said combined image, or a region of 

decreased stiffness being indicated by a region of increased displacement within said 

combined image. 

Step (d) above may optionally be carried out while concurrently delivering an 
interspersed set of pushing pulses to said forcing region to reduce the return of said 
25 target regions from said displaced positions to said initial positions. 

Steps (a) through (e) above may be completed in a total of 50, 25 or 10 
milliseconds or less for each cycle (i.e., each forced region). In some embodiments, a 
cycle of steps (a) through (d) can be completed in 15 milliseconds or less. 

In some embodiments, the pushing pulses are delivered before the first set of 
30 tracking pulses, the initial positions are displaced positions, and the second positions 
are relaxed positions. In another embodiment, the pushing pulses are delivered 
between the first and second set of pulses, the initial positions indicate the relaxed 
positions, and the second positions indicate the displaced positions. 
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Figure 4 illustrates the cyclic repeating of steps (a) through (e) above for 
different forced regions (vertical hatched regions llf, llg, and llj) within the target 
regions (11a through lip) in the axial/azimuthal plane (see Figure 7). The Boxes 
represent the same view as that shown in Figure 6. Arrows represent transition from 
5 one cycle to another (cycles A, B, and C). Note that not all target regions need be 
detected during each cycle, and hence the corresponding transducer elements may be 
active or inactive in various patterns during each cycle. 

As shown in Block 42 of Figure 5, a two-dimensional displacement map can 
then be generated for each cycle A, B, and C of Figure 4. These two-dimensional 

10 displacement maps are then used to generate a single combined image (Block 43) in 
the signal processing device 31 of Figure 3. This combined image can then be 
displayed (Block 44) on the video display device 32 of Figure 3. Of course, the 
single combined image may also be stored in a suitable memory device for future 
reference, printed on a printer, etc. The method may further include generating a B- 

15 mode image of the two dimensional plane in accordance with conventional 
techniques, and then displaying the single combined image superimposed on that B- 
mode image. 

Some embodiments of the invention can be implemented on a Siemens Elegra 
ultrasound scanner, modified to provide control of beam sequences and access to raw 
20 radio frequency data. A Siemens 75L40 transducer may be used as the transducer 
array. 

Some embodiments of the invention may be carried out as follows: 
First, a group of low intensity "tracking lines" that interrogate the tissue 
surrounding the position of interest are fired and stored for tissue initial position 
25 reference. 

Second, a series of one or more focused, high intensity "pushing lines" is fired 
along a single line of flight focused at the position of interest. 

Third, the original group of tracking lines is fired again, in order to determine 
the relative motion caused by the radiation force associated with the pushing lines. 
30 These tracking lines may optionally be interspersed with pushing lines in order to 
reduce or avoid relaxation of the tissue. 

Fourth, each tracking line is divided into sequential axial search regions, and 
the displacements of the tissue within each search region are determined. A number 
of different motion tracking algorithms can be used to determine the relative motion, 
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or displacement, between the initial reference tracking lines and the second set of 
tracking lines fired after radiation force application. Examples include, but are not 
limited to, cross correlation and Sum Absolute Difference (SAD). The a priori 
knowledge of the direction of motion reduces the algorithm implementation time. 
5 Steps 1-4 above are preferably accomplished in 50, 25 or 10 milliseconds or 

less. The results of step 4 are used to generate a two-dimensional displacement map 
of the region of tissue surrounding the position of interest (or force location). 

Fifth, steps 1 through 4 can be repeated, cyclically, for a plurality of force 
locations within a larger two-dimensional imaging plane. The number of forcing 

10 locations and the spatial distribution of the forcing locations may be determined by 
(among other things) the specific transducer, transmit parameters, and the size of the 
region of interest to be interrogated. The same or different sets of elements within the 
transducer array may be used for the tracking pulses with each force location. 

Sixth, each of the two dimensional displacement maps (each of which may be 

15 generated before, during or after subsequent cyclical repeatings of steps 1-4) can be 
combined into a single image (which may or may not be displayed on a video 
monitor, printer or other such display means). Signal processing such as averaging of 
collocated regions, and/or some type of normalization to account for the displacement 
generated in a homogeneous region of tissue, may be employed. 

20 According to some embodiments of the invention, the displacement of the 

tissue may be monitored over time, both while the force is being applied (by 
interspersing the pushing lines and the tracking lines), and after cessation of the high 
intensity pushing lines or pulses. This may be accomplished by firing the group of 
tracking lines repeatedly at the desired time intervals, and evaluating the changes in 

25 the displacement maps over time. 

Figure 6A is an example remote palpation image of a breast lesion phantom. 
This image was generated using multiple forcing locations, separated by 1 millimeter 
(mm) in both the axial and azimuthal dimensions. The image spans a 9 x 1 1 mm 
region, for a total of 99 different pushing locations. Figure 6B provides the 

30 corresponding B-mode image. 

It has been observed that some tissues can exhibit strain-stiffening behavior 
{i.e. glandular tissue, carcinoma) whereas other tissues do not (i.e. fatty breast tissue). 
Therefore, in methods intended to characterize the stiffness of tissue, it is often 
advantageous to pre-compress the tissue. This has the effect of increasing the contrast 
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between the different tissue types (Krouskop et. al., Elastic Moduli of Breast and 
Prostate Tissues Under Compression, Ultrasonic Imaging 20, 260-274 (1998)). This 
will also be the case for the methods described herein. 

For clarity, the interrogation of a two-dimensional plane with multiple pushing 
5 locations (the axial/azimuthal plane - see Figure 7 where transducer array 32 is 
positioned over a target region represented as a cube containing a region of varying 
stiffness 31) has been described. In other embodiments, methods can be carried out in 
a manner that includes the interrogation of a three-dimensional volume. This may be 
accomplished in a variety of ways. According to examples illustrated in Figure 8, 

10 where a transducer array 42 is positioned over a target region represented as a square 
containing a region of varying stiffness 41, and is translated from a first position as 
shown by 42 to a second position shown by 42', one can use the existing planar 
system and translate the transducer in the elevation dimension to sequentially 
interrogate a series of planes comprising a three-dimensional volume. According to 

15 examples illustrated in Figure 9, where transducer array 52 is positioned over a target 
region represented as a square containing a region of varying stiffness 51, one can use 
a two-dimensional transducer array {i.e. one that has several rows of elements), and 
keep the transducer in one location, and steer the beam (represented as lines within 
the cube) to interrogate a three-dimensional sector of the target region. 

20 When using the ultrasound transducer array to either generate the high 

intensity pushing pulses, or the displacement tracking pulses, a set of multiple 
elements may be used to generate each line. The set of elements that is used can 
either comprise all of the elements in the transducer array, or include only a subset of 
the elements. The specific elements that are active for each transmit pulse may be 

25 determined by the desired focal depth, resolution, and depth of field for each line. 
According to some embodiments, the pushing beams can be tightly focused, therefore 
a fairly large number of elements can be used to generate each pushing beam. 

The spatial peak temporal average intensities required to generate detectable 
displacements in tissue vary depending upon the tissue acoustic and mechanical 

30 characteristics. They can be from 10 W/cm 2 to 1000 W/cm 2 , with higher intensities 
being associated with better Signal-to-Noise-Ratios (SNRs). A trade-off may be 
found, however, between increasing intensities and the potential for tissue heating. 
Intensities may be used in the 100 to 400 W/cm 2 range to reduce tissue heating and 
maintain a sufficient degree of intensity. These values may be higher than those used 
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for diagnostic imaging (0.72 W/cm 2 ), and lower than those used for HIFU (High 
Intensity Focused Ultrasound) imaging (1000 W/cm 2 ). Given the short application 
time in a single location, (i.e. 15 milliseconds), the required intensities should not 
pose a significant risk to the patient. 
5 In the method implementation described herein, the high intensity acoustic 

energy can be applied by using a series of multiple, relatively short duty cycle pulses 
(i.e. 40 pulses, each 10 microseconds long, applied over a time period of 10 
milliseconds). The method can also be accomplished by delivering the same amount 
of acoustic energy in a much shorter time period using a single long pulse (i.e. 1 

10 pulse, 0.4 milliseconds long). A sufficient amount of acoustic energy may be 
delivered to the tissue to achieve a given displacement. The amount of energy to 
achieve a given displacement can be accomplished using any number of pulsing 
regimes. One mode of implementation is to use a single, long pulse (i.e. 0.5 
milliseconds), to achieve the initial displacement, and then to intersperse some of the 

15 shorter duty cycle (i.e. 10 microseconds) high intensity pulses with the tracking pulses 
to hold the tissue in its displaced location while tracking. This may reduce the 
amount of time required at each pushing location, and thus reduce the potential for 
tissue heating, while at the same time still achieving the desired tissue displacements. 
The use of a single, long pulse may, however, require additional system 

20 modifications. It may, for example, require the addition of heat sinking capabilities to 
the transducer, as well as modification of a standard power supply to allow the 
generation of longer pulses. 

In some embodiments, tissue heating may be reduced by parallel processing 
techniques, as described in more detail with respect to Figure 10. 

25 The displacement data from each pushing location is combined to form a 

single image. In order to achieve a uniform image, normalization may be useful. 
There are three features may benefit from normalization: 1) attenuation, 2) pushing 
function shape and non-uniformity, and 3) time of acquisition of tracking lines. Each 
of these features may be normalized out of the image, such that a Remote Palpation 

30 image of a homogeneous region of tissue will appear uniform. 

While embodiments according to the invention have been described with 
reference to lesion identification and characterization, some embodiments will also 
find use in any area where information about the varying mechanical properties of 
tissue or any other material subject to ultrasonography may be useful. Some examples 
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include: 1) the characterization of arterial stiffness, which can be indicative of the 
degree of atherosclerotic disease, 2) the assessment of muscle tone, which is of 
importance in determining the course of treatment for female pelvic floor 
complications, and 3) assessing the stiffness of kidneys, which can be indicative of 
5 the viability of kidney transplants. 

Receive mode parallel processing 

As described above, receive mode parallel processing can be used to 
interrogate blood vessels. It will be understood that in some embodiments receive 

10 mode parallel processing can also be used independent of the interrogation of blood 
vessels. In particular, embodiments according to the present invention that can 
receive ultrasound tracking signals from locations in a target region are shown in 
Figure 10. A controller 1035 having a signal analyzer 1031 is operatively connected 
to an ultrasound transducer array 1030. The transducer array 1030 includes 

15 transducer array elements 1030A-1030H that can deliver ultrasound pulses to a tissue 
medium 1032 to interrogate a target region 1036. Reflected ultrasound signals from 
the tissue medium 1032 can be received by the array elements 1030A-1030H. 

As illustrated in Figure 10, an exemplary ultrasound pulse 1034 that is de- 
focused with respect to the target region 1036 propagates from array element 1030B 

20 through the tissue medium 1032 to the target region 1036. Portions of the ultrasound 
pulse 1034 can be reflected towards the transducer array 1030. Reflected ultrasound 
waves can be received from a plurality of directions 1038 by the transducer array 
1030, as shown, by array elements 1030D and 1030H. It should be understood that 
pulses, such as pulse 1034, can be delivered by any of the array elements 1030A- 

25 1030H and that reflected waves can be received by any of the array elements 1030A- 
1030H. 

The signal analyzer 1031 can analyze the various signals received by the array 
elements 1030A-1030H. According to some embodiments of the present invention, 
the signal analyzer 1031 can process the signals received by the array elements 
30 1030A-1030H using receive mode parallel processing techniques. For example, the 
signal analyzer 1031 can sum the signals received by array elements 1030A-1030H 
using a focal delay so that the various signals will be in phase with one another. 
Examples of receive mode parallel processing techniques are shown in U.S. Patent 
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Nos. 5,544,655 to Daigle, 5,718,230 to Chapman et al., and 5,685,308 to Wright et al., 
the disclosures of which are hereby incorporated by reference in their entirety. 

According to some embodiments of the present invention, ultrasound signals 
can be received from locations in the target region, such as is shown in Figure 10. 
5 The received signals can be summed using a focal delay or phase adjustment. For 
example, as shown in Figure 11, a pushing pulse can be delivered from an ultrasound 
transducer array, such as array 1030, to a target region within a medium to displace 
the target region to a displaced position (Block 1100). A tracking pulse can then be 
delivered from the ultrasound transducer array to the target region (Block 1110). A 

10 tracking signal can be received by the transducer array responsive to the tracking 
pulse from locations in the target region (Block 1120). By using receive mode 
parallel processing techniques to analyze the tracking signals from locations in the 
target region, such as by using multiple sums of differently focused delayed signals, 
additional information can be gathered for each pushing pulse so that the frame rate 

1 5 can be increased. In addition, fewer pushing pulses may be used, which can result in 
reduced tissue heating. 

Blood vessel interrogation and/or evaluation 

20 As described above, receive mode parallel processing can be used to 

interrogate blood vessels using ultrasound interrogation techniques. It will be 
understood that blood vessels can be interrogated independently of receive mode 
parallel processing techniques. In particular, further embodiments according to the 
present invention for interrogating a target region 1046 including a blood vessel 1048 

25 are shown in Figure 12. A controller 1044 including a signal analyzer 1041 is 

connected to an ultrasound transducer 1040. The controller 1044 can provide the 
transducer 1040 to deliver a pulse 1045 along a propagation path to a target region 
1046 in a tissue medium 1042. A reflected wave can travel along the propagation 
path to be received by the transducer 1040. The pulses can include ultrasound 

30 "tracking" or "pushing" pulses as described herein, and the reflected wave can be a 
tracking signal. Although a single transducer 1040 is shown in Figure 12, a 
transducer array could also be used and/or receive mode parallel processing 
techniques may be employed. 
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Accordingly, embodiments of the invention maybe useful to determine blood 
vessel health in arterial and veinal vessels, including the popliteal artery and coronary 
artery, as well as heart tissue, and in particular, in the carotid artery of a subject by 
detecting values associated with mechanical properties of a blood vessel or cardiac 
5 tissue. For clarity and ease of presentation, embodiments of the invention may be 
described herein with reference to blood vessels and/or cardiac tissue. However, 
embodiments of the invention may be carried out on blood vessels, cardiac tissue and 
tissue surrounding and/or adjacent a vessel or cardiac tissue. Thus, it is understood 
that a detected mechanical property of a blood vessel or cardiac tissue can include a 

10 detected property of the tissue surrounding the vessel wall or cardiac tissue. 

Advances in imaging and diagnostic techniques, including the techniques 
described herein, may be used to measure certain mechanical properties of a blood 
vessel wall. However, in some cases, it may be difficult to obtain absolute 
measurements of stiffness using such techniques, which may indicate increased 

15 stiffness of a region only relative to surrounding tissue. Furthermore, even if absolute 
measurements of stiffness were easily obtained, it may be difficult to define a clinical 
standard for correlating such a stiffness measurement to vascular health. Each 
individual may have a different baseline healthy level of vessel wall stiffness. An 
amount of stiffness that may be healthy in one individual could be an indication of 

20 disease in another individual, for example, based on age or genetic disposition. 

Numerous variations and implementations of the instant invention will be 
apparent to those skilled in the art. A mechanical property of a vessel wall can be 
detected using methods and systems known to those of skill in the art, including 
detection methods and systems discussed herein and/or disclosed in U.S. Patent No. 

25 6,371,912. As used herein, the term "mechanical property" includes any property or 
characteristic that can be indicative of a mechanical property. For example, stiffness 
is a mechanical property that can be a function of the material that forms an object as 
well as the shape of the object and the tissue surrounding the object. Characteristics 
that can be measured to describe mechanical properties include recovery times, 

30 displacement subsequent to forced movement, time to peak displacement, viscosity, 
and stiffness. In some embodiments, mechanical properties can be measured by 
relative measurements under similar conditions, and therefore, an absolute 
measurement of the mechanical property may be unnecessary to characterize the 
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property. For example, the reaction or vibrations of a material due to the application 
of an ultrasound pulse can indicate mechanical properties of the material. 

With reference to Figure 13, a value associated with a mechanical property of 
a vessel wall and/or cardiac tissue is detected at at least two timing points during the 
5 cardiac cycle of a subject (Block 101). The value of the mechanical property can be a 
measurement of stiffness, viscosity, recovery time constant, displacement, time to 
peak displacement, or a dynamic response to radiation force excitation. The 
mechanical property can be detected at a single point in time or over a period of time, 
typically about 5 ms to about 1 0 ms, For example, an observed response to an 

10 excitation pulse can be observed over a period of time to determine a maximum 

displacement and/or a time to peak displacement. Examples of timing points include 
one point at a relaxed point in the cardiac cycle, and another point at an expanded 
point, or a point during systole and a point during diastole. The values associated 
with the mechanical properties at the two timing points are compared to evaluate the 

15 blood vessel and/or cardiac tissue, such as to determine a measure of vessel and/or 
cardiac health or diagnose a medical condition (Block 103). 

The comparison can be any qualitative or quantitative comparison and can be 
visual or mathematical. The comparison can be carried out by visually observing a 
difference in a mechanical property in an ultrasound image of the vessel wall or 

20 quantitatively detecting a change in the detected mechanical property. The 

comparison may include experimentally determined evaluation standards for health 
and disease that can be determined by a statistical evaluation of a representative 
population. Evaluation standards may be established to account for differences based 
on sex, age, smoking and other factors. The comparison can include normalizing the 

25 evaluation based on vessel distensibility, blood pressure or other physiologic 
characteristics such as sex, age, smoking, and similar factors. The evaluation 
standards can include a threshold value, above or below which is deemed unhealthy or 
a graduated scale for determining a range of relative health. As used herein, a "point" 
can be a one-, two-, or three-dimensional area ranging from an infmitesimally small 

30 point to an area with dimensions on the order of microns, millimeters, centimeters, or 
more. In some embodiments, a "point" can be an isolated region, such as a tissue 
layer in a blood vessel or cardiac area or a focal lesion. 

In other embodiments as shown in Figure 14, a first value associated with a 
mechanical property of a vessel wall and/or cardiac tissue can be detected at a first 
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spatial point within a first layer of a vessel wall (Block 131), and a second value 
associated with a mechanical property at a second spatial point within another tissue 
layer such as a second vessel wall layer (Block 132). The layers may include various 
tissues within the heart or vessel wall layers such as focal lesions, the inner 
5 endothelium layer, the media smooth muscle layer, and the outer adventita connective 
tissue layer. The detected mechanical property is compared at the two spatial points 
to evaluate the blood vessel and/or cardiac tissue (Block 133). The spatial points can 
include areas of vessel wall, including focal lesions, and/or cardiac tissue. Preferably, 
the two tissue layers are the media and adventita layers. 

10 Without wishing to be bound by theory, it is currently believed that the relative 

mechanical properties of diseased vessel tissue during various points in the cardiac 
cycle or at various points within different vessel wall layers may be used to measure 
blood vessel health. Mechanical properties at a single spatial or temporal point may 
provide some indication of blood vessel health, however, the healthy baseline of a 

1 5 mechanical property such as stiffness, recovery time constants, viscosity, or 

displacement measurements may vary based on age, gender, ethnicity, location in the 
body, or even randomly between individuals. Thus, comparing relative mechanical 
properties during various points in the cardiac cycle as the vessels expand and 
contract or between vessel wall layers may provide a measurement of vascular health. 

20 Such relative measurements in a subject may account for varying individual 

characteristics. It is currently believed that a greater relative change in a mechanical 
property such as stiffness indicates lesser arterial health over smaller relative changes. 

As discussed briefly above, examples of detected mechanical properties 
include a dynamic response to radiation force excitation, a measurement of stiffness, a 

25 measurement of a recovery time constant, a measurement of viscosity, a measurement 
of time to peak displacement and/or a measurement of displacement. Mechanical 
properties can be derived measurements based on observed responses of tissue to an 
excitation force. Various parameters can be selected from the observed responses as a 
detected mechanical property without necessitating an absolute measurement or direct 

30 measurement of a mechanical property. For example, mechanical properties can be 
determined from an evaluation of the dynamic response of tissue to radiation 
excitation forces. The measurement can be made at a single point in time or over a 
period of time (typically about 5 ms to about 10 ms), for example, to determine the 
time to peak displacement. The mechanical property of the vessel wall can be 
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detected at one point during systole and at another point during diastole of the cardiac 
cycle. It is believed that detection of the mechanical property during systole and 
diastole can maximize the difference between the mechanical property at the two 
points, thus increasing the sensitivity of a measurement, because measured difference 
5 in the detected mechanical property is between a relaxed blood vessel wall and an 
expanded vessel wall. 

With reference to Figure 15, further embodiments include detecting a value of 
a mechanical property of a vessel wall and/or cardiac tissue at two spatial points 
within two layers, such as two cardiac or vessel wall layers (Block 141). The value of 

1 0 the mechanical property can be detected at the two spatial points at at least two timing 
points during the cardiac cycle (Block 143). The detected values of the mechanical 
property at the two spatial points and at the two timing points may then be compared 
to evaluate the blood vessel and/or cardiac tissue (Block 145). In some embodiments, 
one mechanical property measurement is made during systole of the cardiac cycle and 

1 5 another mechanical property measurement is made during diastole of the cardiac 
cycle. 

The mechanical property can be detected by methods and systems known to 
those of skill in the art capable of detecting mechanical properties of vessels as 
described herein, including remote palpation, ARFI, and receive mode parallel 

20 processing ultrasound techniques. For example, systems may be provided for 
measuring the mechanical property multiple times or at multiple spatial points using 
ultrasound measurements of the blood vessel wall and comparing those measurements 
to evaluate blood vessels. 

For example, with reference to Figure 12, beam sequences, which can be 

25 delivered by transducer 1040, can be repeated throughout the cardiac cycle. The 
controller 1044 can be configured to gate the timing of the radiation force 
interrogations delivered by the transducer 1040 with an electrocardiogram. Other 
ways of timing the radiation force interrogations or determining the timing of the 
collected data can be used. For example, displacement maps can be generated over a 

30 period of time, and the timing of the map with respect to cardiac cycle can be 
determined from the map itself. The radiation force interrogations can also be gated 
with a blood pressure measurement to determine the timing of the cardiac cycle. Data 
can be collected by the signal analyzer 1041 at specific, predetermined points during 
the cardiac cycle. Multiple locations along the vessel wall and/or cardiac tissue, 
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including the surrounding tissue, can be evaluated, either using simultaneous 
measurements or by moving the transducer to evaluate different locations. Data can 
also be collected by the signal analyzer 1041 throughout a series of closely spaced or 
continuous measurement during the cardiac cycle. Thus, the entire cardiac cycle can 
5 be observed. 

Displacement magnitude is generally inversely proportional to local tissue 
stiffness and, in some embodiments, may be on the order of ten microns. The 
transient shear wave velocity and shear wave attenuation are also related to the tissue 
mechanical properties. The volume of tissue to which radiation force is applied can 

10 be determined by the focal characteristics of the transmitting transducer. Comparison 
of detected mechanical properties can be accomplished in various ways. For example, 
the ratio of matched parameters such as maximum displacement or recovery time 
constant can be compared. Such comparisons may be made at various points in the 
cardiac cycle, such as at systole and diastole, and/or with respect to measurements at 

1 5 different blood vessel wall layers or over a period of time. 

Figure 16 illustrates operation of the system of Figure 12 in further detail. As 
seen in Figure 16, for each of the spatial or timing points (Block 111) a set of one or 
more tracking pulses are delivered from a plurality of transducer elements in an 
ultrasound transducer array, such as transducer array 1040 in Figure 12, to a plurality 

20 of target regions in a two-dimensional plane intersecting the vessel wall to detect 
initial positions for target regions (Block 113). A first set of one or more pushing 
pulses is delivered from the transducer array 1040 to a forcing region among the 
target regions to displace the target regions (Block 115). A second set of one or more 
tracking pulses is delivered from the transducer array 1040 to the target regions to 

25 detect displaced positions (Block 117). The step of Block 117 can be repeated for a 
single set of pushing pulses over a period of time, such as between about 5 ms and 
about 10 ms. The repetition of the tracking pulses in the step of Block 117 can be 
used to observe a reaction over a period of time, for example, to determine 
displacement over time and/or the time to peak displacement for an interrogation 

30 pulse. In addition or alternatively, the steps of Block 113, Block 115, and Block 117 
may be repeated at the same point. The steps of Block 113, Block 115, and Block 
117 may be repeated one or more times for the next point (Block 119). The points at 
which the mechanical property is measured may be timing points at different points 
during the cardiac cycle and/or spatial points at different positions within blood vessel 
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walls. In some embodiments, the steps of Block 113, Block 115, and Block 117 can 

be repeated in a series of cycles, with the pushing pulses being delivered to a different 
forcing region and diferent target regions during each of the cycles. 

For each of the points, a two-dimensional displacement maps from each of the 
5 initial positions and displaced positions can be generated to produce a plurality of 
two-dimensional displacement maps (Block 121). The plurality of two-dimensional 
displacement maps can be made over time to show characteristics such as 
displacement over time. The plurality of two-dimensional displacement maps can be 
combined into a single image for each point, with a region of increased stiffness being 

10 indicated by a region of decreased displacement (Block 123). The images can be 
combined into a single parametric image, for example, isolating one or more 
evaluation criteria. Two or more spatial and/or timing points can be compared to 
evaluate the blood vessel and/or cardiac tissue (Block 125). Images may also be 
made in one or three dimensions. 

15 The operations shown in Figure 16 can be used to detect mechanical properties 

for both timing and spatial points at various times during the cardiac cycle and/or 
within two or more tissue layers, such as cardiac tissue layers and/or vessel wall 
layers. 

Other imaging techniques and non-imaging data analysis techniques known to 
20 those of skill in the art may be used to compare spatial and timing points at which a 
mechanical property is detected. 

Figures 17 and 18 are images produced using ultrasound techniques described 
in U.S. Patent No. 6,371,912. Figure 17 is a displacement image obtained by 
operations shown in Figure 16 of an in vivo carotid obtained during diastole, which 
25 shows a delineation of the posterior vessel wall. Figure 18 is a magnified illustration 
of the image in Figure 17. Figure 19 illustrates the displacement value (plotted on 
the horizontal axis) as a function of position from the data used to produce Figures 17 
and 18. The images in Figures 17, 18, and 19 were generated using displacement 
data 0.1 msec after radiation force application. 
30 While the present invention may be used to detect regions of varying stiffness 

{i.e., increased or decreased elastic modulus) in any type of medium, the medium is in 
a preferred embodiment a biological tissue such as breast tissue in a living subject, 
where the regions of varying stiffness to be detected are tumors. 
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The invention may be carried out on human subjects for diagnostic or 
prognostic purposes, and may be carried out on animal subjects such as dogs and cats 
for veterinary purposes. 

Circuits and devices described herein may be implemented as general or 
5 special purpose hardware devices, as software running on or programmed into general 
or special purpose hardware devices, or as combinations thereof 

Numerous variations and implementations of the instant invention will be 
apparent to those skilled in the art. Conventional ultrasound apparatus is known, and 
is described in, for example, U.S. Patent No. 5,487,387 to Trahey et ah; U.S. Patent 

10 No. 5,810,731 to Sarvazyan and Rudenko; U.S. Patent No. 5,921,928 to Greenleaf et 
al.; M. Fatemi and J. Greenleaf, Ultrasound-stimulated vibro-acoustic spectrography, 
Science, 280:82-85, (1998); K. Nightingale, Ultrasonic Generation and Detection of 
Acoustic Streaming to Differentiate Between Fluid-Filled and Solid Lesions in the 
Breast, Ph.D. thesis, Duke University, 1997; K. Nightingale, R. Nightingale, T. Hall, 

15 and G. Trahey, The use of radiation force induced tissue displacements to image 
stiffness: a feasibility study, 23 rd International Symposium on Ultrasonic Imaging and 
Tissue Characterization, May 27-29, 1998; K.R. Nightingale, P. J. Kornguth, S.M. 
Breit, S.N. Liu, and G.E. Trahey, Utilization of acoustic streaming to classify breast 
lesions in vivo, In Proceedings of the 1997 IEEE Ultrasonics Symposium, pages 

20 1419-1422, 1997; K.R. Nightingale, R.W. Nightingale, M.L. Palmeri, and G.E. 
Trahey, Finite element analysis of radiation force induced tissue motion with 
experimental validation, In Proceedings of the 1999 IEEE Ultrasonics Symposium, 
page in press, 1999; A. Sarvazyan, O. Rudenko, S. Swanson, J. Fowlkes, and S. 
Emelianov, Shear wave elasticity imaging: A new ultrasonic technology of medical 

25 diagnostics, Ultrasound Med. Biol 24:9 1419-1435 (1998); T. Sugimoto, S. Ueha, and 
K. Itoh, Tissue hardness measurement using the radiation force of focused ultrasound, 
In Proceedings of the 1990 Ultrasonics Symposium, pages 1377-1380, 1990; and W. 
Walker, Internal deformation of a uniform elastic solid by acoustic radiation force, J. 
Acoust. Soc. Am. y 105:4 2508-2518 (1999). The disclosures of these references are to 

30 be incorporated herein by reference in their entirety for their teaching of various 
elements and features that may be used to implement and carry out the invention 
described herein. 
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The foregoing is illustrative of the present invention, and is not to be construed 
as limiting thereof. The invention is described by the following claims, with 
equivalents of the claims to be included therein. 
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